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The Gaver-Stehfest algorithm for numerical inversion of Laplace transform was developed in the 
late 1960s. Due to its simplicity and good performance it is becoming increasingly more popular 
in such diverse areas as Geophysics, Operations Research and Economics, Financial and Actuarial 
Mathematics, Computational Physics and Chemistry. Despite the large number of applications and 

numerical studies, this method has never been rigorously investigated. In particular, it is not known 
whether the Gaver-Stehfest approximations converge and what is the rate of convergence. In this 
paper we answer the first of these two questions: We prove that the Gaver-Stehfest approximations 



converge for functions of bounded variation and functions satisfying an analogue of Dini criterion. 
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1 Introduction and main results 



In this paper we are concerned with the classical problem of numerical inversion of Laplace transform. 
More specifically, assume that / : (0, oo) i— M is a locally integrable function, such that its Laplace 
transform 



oo 

F{z) ■= j e-^V(x)dx 



is finite for all 2; > 0. The problem consists in recovering the original function f{x) given that we know 
F{z). This problem has numerous applications, and it has attracted a lot of attention from researchers 
over the last fifty years (see [3] for an up-to-date exposition of this area). 

More specifically, we are interested in the Gaver-Stehfest algorithm, which aims to approximate f{x) 
by a sequence of functions 



2n 



fn{x) := ln{2)x-^ J] ak{n)F {-n ln(2)a;-^) , n>l, x>0, (1) 



k=l 

where the coefficients are defined as follows: 

.n+k min(fc,n) 



ak(n 



In order to demonstrate the intuition behind these rather complicated formulas, let us explain where 
they come from. In 1966 Gaver [8] has introduced the following sequence of approximations 

Ux) := ln(2)x-^^^^^ g Q {-iyF{{k + i) H2)x~'), k>l,x>0. (2) 

Applying the binomial theorem, the above expression can be rewritten in an equivalent integral form 

00 

~ f f xu \ 

fk{x) = / Pk{u)f —- du, (3) 







ln(2) 



where 



Pk{u) ^= fe,(fc!.^i), (l-0'e-^ k>l,u>0. 

The function Pk{u) is positive and its integral over [0, 00) is equal to one, therefore we can think of it 
as the density function of a positive random variable Uk- One can check by direct calculations that 
E[[/fe] = ln(2) + 0{k~^) and Vai{Uk) = 0{k~^) as k ^ +00. This shows that the random variables 
{Uk}k>i converge in distribution to ln(2), therefore for any continuous function f{x) and any x > we 
have 



fkix) = E 



f 



xUi 



ln(2))j 



— i- f{x), A; — !■ 00. 
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Gaver [8] also proves that 



AW = /W + ^ + ^ + ... + ^ + o(r-^). fc^+oo. (4) 

under the assumptions / G C'^"^(M.~^) and f^^\x) G Loo(K^) for all j = 0, 1, . . . , 2m. Formula (4) shows 
that fk{x) converge to /(x) rather slowly, but it also suggests using a convergence acceleration technique. 
This is exactly what was done by Stehfest [18] in 1970. He defined a new approximation 

n 

fn{x) = ^Ck{n)fk{x), (5) 
fc=l 

where the coefficients Ck{n) are chosen so that the asymptotic terms bj{x)k~^ , j = 1, 2, . . . , — 1 in (4) 
are eliminated: 

y Ckin)k-^ = l'^' (6) 
tl 10' J = l,2,...,n-1. ^ ' 

It turns out that the above conditions are satisfied by the sequence 

c,{n) :=(-l)"+^ ^" n> 1,1 <k<n. 

k\[n — k)\ 

Combining the above formula with (2) and (5) gives us the Gaver-Stehfest approximations in the form 
(!)• 

The Gaver-Stehfest algorithm has a number of desirable properties: (i) the approximations fn{x) are 
linear in values of F{z); (ii) it requires the values of F{z) for real z only does not need complex numbers 
at all; (iii) the coefficients ak{n) can be easily computed; (iv) the Gaver-Stehfest approximations are 
exact for constant functions, that is, if f{x) = c then fn{x) = c for all n > 1. This algorithm was studied 
in [6, 10, 15, 19], where it was demonstrated numerically that fn{x) converge very quickly to f{x) for 
many examples of initial functions f{x) (provided that f{x) is non-oscillating). Another universally 
accepted fact is that this algorithm requires high-precision arithmetic for its implementation (which is 
rather obvious, since the coefficients ak{n) are growing very rapidly and have alternating signs). 

Over the last forty years the Gaver-Stehfest algorithm has been applied to solving various problems in 
Geophysics [13], Probability [1, 14], Actuarial Mathematics [2] and Mathematical Finance [17], Chemistry 
[16] and Economics [11]. This is just a small sample, and a quick search on the Internet produces hundreds 
of papers with references to the Gaver-Stehfest algorithm. 

Despite all this popularity, there has not yet been a single rigorous investigation of this algorithm. In 
particular, it is not known what are the sufficient conditions on / that will ensure convergence of fn{x), 
and what would be the rate of convergence. Stehfest [18] writes that "theoretically fn{x) become the 
more accurate the greater n", and Proposition 8.2 in Abate and Whitt [1] states that for any A; > we 
have fn{x) — f{x) = o{n~^) as n — )■ +cxd, however both of these statements are not very precise and lack 
rigorous proof. The authors seem to assume that the derivation of the Gaver-Stehfest approximations 
via formulas (2), (4), (5) and (6) automatically gives the proof of convergence. This is not the case. First 
of all, validity of the asymptotic expansion (4) for all m > 1 would require a very restrictive assumption 
on the function / (it has to be an infinitely differentiable function plus some additional assumptions). 
The second issue is that even if the asymptotic expansion (4) is valid for all m > 1, it is still not clear 
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that fn{x) — )■ f{x), since the coefficients Cfc(n) grow very rapidly and have alternating signs. Therefore, 
it is not at all obvious that fn{x) will converge to f{x), given (5), (6) and the fact that fn{x) — )■ f{x). 

In this paper we present the first rigorous investigation of the Gaver-Stehfest algorithm. We derive 
two sufficient conditions on the function /, which ensure convergence of fn{x)- Our main results are 
presented in the next theorem. 

Theorem 1. Assume that f : (0, oo) i— )■ M zs a locally integrable function, such that the Laplace transform 
F{z) = e~^^ f{x)dx exists for all z > 0, and that fn{x) are defined by (1). 

(i) The convergence of fn{x) depends only on the values of the function f in the neighborhood of x. 

(a) Assume that for some c G M and some e G (0, 1/4) 

j |/(-xlog2(^ + v)) + f{-x\og^{l - v)) - 2c\v-^dv < oo. (7) 



Then fn{x) c as n ^ +oo. 

(Hi) Assume that the function f has bounded variation in the neighborhood of x. Then 
fnix) ^ (fix + 0) + fix - 0))/2 asn^ +oo. 

Our sufficient conditions are very similar to the corresponding results from the theory of convergence 
of Fourier series. Item (i) has its counterpart in [12, Theorem 4.1.1], item (ii) should be compared with 
Dini criterion [12, Theorem 4.1.3(i)] 

e 

j \f{x + v) + fix -v)- 2c\v'^dv < oo, (8) 



and item (iii) is exactly the same as Jordan criterion [12, Theorem 4.1.3(ii)]. Theorem (l)(ii) also provides 
the following useful corollary: 

Corollary 1. Under the assumptions of Theorem 1, if 

fix + v)-fix) = Oi\vn 

for some a > and all v in some neighborhood of x, then fn{x) — )■ f{x) as n ^ +oo. 

We do not address the second important problem related to the rate of convergence of fn{x)- We 
hope that the methods developed in this paper can be useful for solving this problem, and we leave it 
for future work. 

2 Proof of Theorem 1 

Let us review some properties of the Lambert W-function (see [5]), which will be used in the proof of The- 
orem 1. We will be interested in the principal branch of the Lambert W-function, which will be denoted 
by W{z)): It is an analytic function in the neighborhood of 2; = and it satisfies W{z) exp(W{z)) = z. 
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Figure 1: Domain (b) and range (a) of the principal branch of the Lambert W-function. 



Let us investigate the range and domain of W{z) in more detaiL The function w = x + iy ^ we^ 
takes real values only if ?/ = or x = —ycot{y) (see Section 4 in [5]). In particular, the open set 

A := {w = X + iy, x > —ycot{y), — tt < y < n} 

is mapped hj z = we^ onto the cut complex plane C\(— oo,— 1/e], see Figure 1. The function w i-)- we'^ 
is one-to-one on the set A, and W{z) is defined as the inverse of this function. Therefore, W{z) is 
analytic in C \ {—oo, —1/e] and maps this set onto A. It is known that W{z) has the following explicit 
Taylor series at ^ = (see formula (3.1) in [5]) 



W{z) = J2{-nr-'-j, \z\<l/e, (9) 



n>l 



and has a branching singularity at 2; = —1/e 

TT./ N V- n . p2 11 3 43 4 769 5 
W(z)=y u„p" = -1 + n - — H p^H P + (10) 



n>Q 



where p = p{z) := 1^2(1 + ez) and the series converges for \p\ < y/2 (see formula (4.22) in [5]). 

We will extend function W{z) to the interval G M, z < —1/e} so that W{z) is continuous in 
the upper half-plane lm{z) > 0. Thus W{z) maps the interval {2; G M, 2; < — 1/e} onto the curve 
X = —ycot{y), y G (0,7r) (which is the upper half of the black curve on Figure 1(a)). The function W{z) 
is smooth on {2; G M, z < —1/e}, and it is clear that |Vr(2;)| and lm(W{z)) are decreasing function of 
z. Assuming that z is real and z < —1/e, the equation we"^ = z has infinitely many solutions, which 
correspond to different branches of Lambert W-function. However, only two of these solutions (given by 
w = W{z) and w = W{z)) satisfy |Im(iy)| < tt. These solutions necessarily lie on the black curve on 
Figure 1(a), and they will play a very important role in our investigation. 

Our main tool used in the proof Theorem 1 is the following sequence of polynomials 
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Recall that {a)k '■= (i{a + 1) ... (a + A; — 1) denotes Pochhammer symbol. Combining equations (3) and 
(5) we obtain an integral representation 



fn{x) 



oo 

I {Ae--{1 - e-)) / du, (12) 



which explains why these polynomials are important for Gaver-Stehfest approximations. Our plan for 
proving Theorem 1 is to establish uniform asymptotics for g„(f), < v < 1, as n — )■ oo and then use 
(12) to study convergence of fn{x). 
Let us also define 

n>l ^^-^ 

Using Stirling's asymptotic formula for the Gamma function one can check that the above series converges 
for \z\ < 1/e. 

Proposition 1. For < < 1 and \t\ < l/(2e) 

G(t;te*) = J]g„(t;)(-l)"r. (14) 

n>l 

Proof. Using (11), the trivial estimates {^)k < {l)k = k\ and k^^^ < n"^^ and the binomial theorem, we 
obtain an upper bound 

£ E (^r^I^ < all < „ < 1. 

which implies that the series in the left-hand side of (14) converges absolutely for \t\ < l/(2e). We 
combine (11) and (14), interchange the order of summation and obtain 



n>l n>l 



J2^^i-lfk'^\vter = G{vte% 
>i ^ 



□ 



It is well known that the asymptotic behavior of the coefficients of the power series is closely related 
with the asymptotic behavior of the generating function at its dominant singularities. Our next goal 
is to obtain analytic continuation of the function G{z) and to study its asymptotic expansion at the 
dominant singularity. 



Proposition 2. There exists a function g{z), which is analytic in C\' 
half-plane lm{z) > 0, such that 



-oo, —1/e] and continuous in the 



G{z) 



V^TT 



(l + e^)"' + ^Ml + e^) 



+ g{z). (15) 
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Proof. Let us denote H{z) := — (^z-^Y^i^)- From (9) we find that 

^ , „n+l 

H{z) = Y.—{-irz-, \z\<l/e. 

n>l 

Using the following identity (see formula 3.621.3 in [9]) 



(sin(t))""dt 
■n J n\ 


and formula (13) we obtain the integral representation 

2 

G{z) = ^ j H(zsm(tf)dt, \z\ < 1/e. (16) 


The function W{z) is analytic in the cut plane C \ (—00, —1/e], therefore H{z) is analytic in the same 
domain. Applying Leibniz rule to (16) we see that G{z) is also analytic in C \ (—00, —1/e]. 
From (10) and the above definition of H{z) we obtain 

H(z) = y Cnpizr = p(z)-^ - —p(z)-^ _ A _ + (17) 
K J nF\ ) F\ ) 24^^ ^ 135 1152 ' ^ ' 

n>-3 



where p{z) := a/2(1 + ez) and the above series converges for \p{z)\ < \f2. We subtract the dominant 
asymptotic terms from -ff(;z) and define the two functions 

TT 

2 

M^):=/f(z)-p(z)-3 + ^p(^)-i = J]c„p(z)", 9i{z):='^- f h{zsmitf)dt. (18) 

It is clear that h{z) is analytic in C \ (—00, —1/e] and continuous in the upper half-plane Im(z) > 0, 
therefore gi{z) satisfies the same properties. 

Next, we use formula 3.681.1 in [9] and formula (14) on page 110 in [7] and conclude that for all 
z G C\ (-00,-l/e] 



2 

2 /* 3 2 1 

- / (l + e;2sin(t)2)-idt = 2Fi{l '^■,1] -ez) = —— - —ln{l + ez) + g2{z) 

TT J Ti[i + ez) Zn 



(19) 



2 

- [{l + ezsm{t)^)-'^dt = ^F,m;l;-ez) = --\n{l + ez)+gs{z), (20) 

TT J TT 


where 2-^1 (c^i b; c; z) denotes the hypergeometric function and both functions (72 (-2) and (73 (-2) are analytic 
in C \ (—00, —1/e] and continuous in the half-plane lm{z) > 0. Combining (16), (18), (19) and (20) we 
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see that 



G{z) 



2 2 2 

2 f 11 2 f 2 f 

H{z?>m{tf)dt= - / p{zsm{tf)~^dt x - / p{z sm{tf)-^ &t + - / h{zsm{tf)dt 

TT J 24 vr J IT J 



+ 



V27t{1 + ez) 24v^7r 
which is equivalent to (15). 



ln(l + ez) + [giiz) + g^iz) + g^{z)] , 



□ 



Lemma 1. For any e G (0, 1) there exist constants b = 6(e) > 1 and C = C{e) > such that 
\qniv)\ < Cb-'^v for all v e[0,l-e]. 

Proof. The function te* maps the unit circle |t| < 1 into some open domain which is a subset of C \ 
(— oo, — 1/e] (see Figure 1). Using this fact combined with the results of Proposition 2, we see that the 
function t G{te^) is analytic for |t| < 1, therefore there exists R > 1 such that t h-)> G{{1 — e)te*) is 
analytic for \t\ < R. It is clear that for all v E [0, 1 — e] the function t i— )■ G'{vte*)te* is also analytic for 
|t| < R. 

We differentiate both sides of (14) with respect to v and find 



{-ly 



d" 



dt 



- {G\vte')te') 



We set b = {1 + R)/2 and use Cauchy estimates (see [4, Theorem 2.14]) 

d'^ 



dr 



{G'{vte')te') 



t=o 



where := max{|G'(fte*)te*| : \t\ = b}. Applying the maximum modulus principle shows that for all 
G [0, 1 - e] 

< max{|G"(t;te*)| : |t| = 6, < i; < 1 - e} x max{|te*| : \t\ = b} 
= max{|G'((l - e)te*)| : \t\ = b} x be^ = C. 

So far we have established that there exist 6 > 1 and C > such that |q',^(f)| < C6~" for all v E [0, 1 — e]. 
Since gn(0) = 0, we conclude that |(3'n('?^)| < Jq |Q'n(^)|d^ < Cb~"-v for all f G [0, 1 — e]. □ 



From now on, we will use the notation w = w{v) := W{—l/{ev)) where v G (0, 1] and W is principal 
branch of the Lambert W-function. In other words, z = w is the unique solution of the equation 
1 + vze^~^^ = in the strip < Im(2;) < tt (see the discussion on page 5, preceding formula (11)). Note 
that w{l) = —1 and both lm{w{v)) and |w(f)| are decreasing functions of f G (0, 1] (see Figure 1). 



Proposition 3. As n +oo we have uniformly for all v E [1/2, 1) 



Qn[V) 



-l)"^Re 

TT 



W~ 



(l + w) 



+ 0(u;-"). 



(21) 
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Proof. As was noted in the proof of Lemma 1, for any v e [0, 1] the function z h- )■ G{vze^) is analytic in 
the circle \z\ < 1. Therefore, using Cauchy integral formula and (14) we find that for all v G [1/2, 1) 



2m 



G{vze^)z 



-ra-l 



dz 



(22) 



Ul/2 



where the contour of integration Ur denotes the circle of radius r and center at 2; = 0, winding counter- 
clockwise around the origin. 

Let g{z) be the function defined by (15). From (15) and (22) we find 



qn[v) 



where we have defined 



Qn,l{v) 



Qn,2[V} 



V27r 



1 

27ri 

1 
27ri 

1 
27ri 



/2 



24v^7r 

1 + vze^~^^ 
ln(l + vze 



-<ln,2{v) + qn,z{v) 



(23) 



g{yze'^)z 



-n-l 



dz. 



Ul/2 



Our first goal is to prove that for all v G [1/2, 1) we have 



(ln,l[V) 



2Re 



'l+w) 



+ 0{w "), n +00, 



(24) 



and that the implied constant in 0{w ") does not depend on v. One can check numerically that 
|w(l/2)| = \W{-2/e)\ ^ 1.2508 < 2. In particular, for all v G [1/2,1) we have 1 < \w{v)\ < 2 and 
the equation 1 + vze^^^ = has exactly two solutions in the circle U2, given by z = w and z = w, 
and no other solutions in the circle 11^^. This shows that for v G [1/2,1) the meromorphic function 



z i—^ z' 



-n-l 



/(I + vze ~^^) has two simple poles at z = w and z = w and no other singularities in the circle 
Ut^, while if f = 1 we have a double pole at z = w = w = —1. We assume that v G [1/2, 1), compute the 
residues of z~^^^/{l + vze^^^) at z = w and z = w and applying Cauchy's residue theorem we obtain 



1 
27ri 



1 + vze^~^^ 



dz = 2Re 



w " 






.(1 + ^). 




1 + f ze^+^ 



dz 



Ul/2 



Us 



The integral in the right-hand side of the above formula is estimated as follows 



1 
27ri 



1 + vze^~^^ 



dz 



< 



3-n 



0{w-n, 



where Ci := min{|l + vze^^^\ : 1/2 < v < l,\z\ = 3}. Note that Ci is strictly positive, since we have 
shown above that for v G [1/2, 1] the function z 1 + vze^^^ has no roots in the domain 2 < < vr. 
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Our second goal is to prove that for all f G [1/2, 1) we have 

qn,2{v) = 0(w""), n +00. 
We isolate the dominant singularities at z = w and z = w and obtain 



(25) 



1 

27ri 



J z-"-Mn(l - z/w)d^ + ^ j z-^'-Hnil- z/w)dz 



Ul/2 



Ul/2 



+ 



27ri 



-n-l 



In 



Ul/2 



1 + vze^^^ 



(1 — z/w){l — z/w) 



dz. 



The first and the second integrals can be evaluated explicitly to (— n^^w~") and (— n^^w^") respectively, 
thus they contribute 0(w~"). The integrand in the third integral is analytic in the circle \z\ < it, and 
using the same technique as was used above for estimating we find that for all v G [1/2, 1) the 

third integral is bounded from above by 0(3~") = 0{w~"'). 

Finally, for any v G [1/2, 1] the function z 1— )■ g{vze^) is analytic in the circle U\yj\ and continuous on 
the boundary of this circle (note that g{z) is continuous at z = — 1/e (see Proposition 2). Therefore 



2m 



g{vze'')z '"■ dz 



27ri 



g{vze'')z '"■ dz 



Ul/2 



I ^ I 



and we obtain the estimate 



\<ln,z{v) 



1 

27ri 



g{vze'')z " dz 



Uu 



< |m;|-"Co 



where C2 '■= ma:x{\g{vze^)\ : 1/2 < f < 1,2; = jwle^'^^*, < t < vr}. Since w{v) is continuous for 
V G [1/2, 1] and g{z) is continuous in the upper half-plane lm{z) > 0, we see that C2 is finite and we 
obtain 



= 0{w "), n +00. 



(26) 



for all V G [1/2,1). 

Combining (23), (24), (25) and (26) gives us (21). 



□ 



Remark 1. With some extra work one can improve the results of Propositions 2 and 3 in the following 

way: 

(i) The function G{z) admits an asymptotic expansion at 2; = — 1/e of the form 



G{z) 



V^TT 



(1 + ez)-' + ^(1 + ez)"(a„ ln(l + ez) + 6„) 



n>0 



where oq = 5/24, ai = 25/1152 and {an}n>2 are computable rational numbers. 
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(ii) As ri — )■ +00 we have uniformly in v G [1/2, 1) 

5 



= (-l)"^Re 

71 



W 



25 



-w 



[1 + w)w' 



l + w) 24n ' 1152^2 
where = W{—l/{ev)). The above formula remains valid as f — )■ 1~, in which case we obtain 



gn(i) 



V2 



7T 



1 



n 



3 2An 



+ 0{n-^). 



For t> G [0, 1/2) we define 

iiy) = w{l - Av^) = W^(-l/(e(l - Av'^))), a{v) = Im [^(f )] . 

These two function will play an important role in the proof of Theorem 1. We summarize some of their 
properties in the next lemma. 

Lemma 2. 

(ii) The function \^{v)\ is smooth and strictly increasing for v G [0, 1/2). It is analytic in the neighborhood 
of V = and satisfies 



8 14 /2 

^(v) = -1 + 2V2W + -v^ + -^iv^ + 0{v^), V ^ 0+. 

o 9 



(27) 



(ii) The function a{v) is smooth and strictly increasing for v G [0, 1/2). It is analytic in the neighborhood 
of V = and satisfies 
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(Hi) For any function h G Li{0, 1/2) we have 

1 
2 

lim / h{v)^{vydv = 0. 



(2^ 



Proof. Items (i) and (ii) follow easily from properties of Lambert W-function. In particular the series 
expansion (27) follows from (10), while (28) is a simple corollary of (27). Item (iii) is obvious, given that 
1^(0)1 = 1 and \^{v)\ is a strictly increasing function. □ 



Proposition 4. Under the assumptions of Theorem 1, fn{x) c as n 00 if and only if for any 
e G (0, 1/4) 



e 

lim / \^{v) 

n— >oo / 



_„sin(na(f)) 



aiv) 



[f {-xlog, {l+v))+f {-xlog, il-v))- 2c] dv = 0. 



(29) 
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Proof. Assume that e G (0, 1/4) and define uq := — ln(^ + e) and ui := — ln(^ — e). Note that < mq < 
ln(2) < ui < oo. The function u i— )■ 4e~"(l — e~") is strictly increasing for < m < ln(2) and strictly 
decreasing for ln(2) < u, and its value at m = ln(2) is one. We conclude that there exists 6 = 5(e) > 0, 
such that for all u G [0, uq] U [ui, oo) we have < 4e~'"(l — e""") < 1 — 5. Using Lemma 1 we find that 
there exist 6 > 1 and C > such that 

\qn (4e-"(l - e"")) | < 4Cr"e-"(l - e"") < 4Cr"e-" 

for all u G [OjMo] U [ui, oo). This fact and our assumptions on the function / guarantee that 



[0,«o]U[ui,oo) 

Thus we have established the following result: 

Ml 

fn{x 



J qn (4e-"(l - e-")) / (^^^ du + o(l), n ^ +oo. 



UO 



As we have pointed out in the introduction, Gaver-Stehfest approximations are exact for constant func- 
tions: if f{x) = C then /„(z) = C for all n > 1. This result and the above formula imply 



fn{x) 



/.«(4e-"(l-0) [/(^) 



dw + ofl), n — )■ +00. 



(30) 



Our next goal is to simplify the integral in (30). We separate this integral into two parts: the integral 
over [Mo,ln(2)] and the integral over [ln(2),Mi]. For the first {resp. second} integral we change the 
variable of integration u = — ln(| + v) {resp. u = — ln(^ — v)} and after combining the two parts we 
obtain 



f^{x)-c= / qn{l-^v^) 



f{-x\0g2il+v)) - C f{~x\0g2il-v)) - C 



+ 



dv + o(l), n +00 (31) 



For our next step, we will need the following result: 



,,.(i-4,.=) = ^kwr"5'=^+oKw 

TT a[v) 



(32) 



uniformly for all v G (0, 1/4]. Note that the Lambert W-function satisfies W{z) = z exp{—W{z)), thus 
for z < — 1/e we have ajLg{W{z)) = vr — Im[14^(2;)], which implies arg(^(t')) = vr — a{v) This gives us 

^(t;)-" = |^(^)|-«e-*'^('^-°(^)) = (-l)"|^(t;)|-" (cos(na(t;)) + i x sm{na{v))) . 

Combining the above formula with the identity 

1 1 + ^ 1 + Re(0 _ . Im(0 



1 + e |l + ep + 
12 



we obtain 



Re 



l + Re(^(f)) a(v 
cos[na[v))— — h sm(na(t')) 



Formulas (27) and (28) show that uniformly for all v G (0, 1/4] 

a(v) 1 



ii^^o(i). 



thus 



Re 



|l + e(^^) 
e(^)-" 



0(1) 



a[v) 



which together with (21) implies (32). 
Next, we define 



g{x,v,c) := / (-xloga (I + w)) + / (-xlog2 (I - t;)) - 2c 
g{x,v,c) : = 



/ (-a; logs +^) ) - c / (-xlog2 (I -t;)) -c 

~r 1 



It is easy to check that 



+ v 



+ 



2g{x,v,c) + 2vg{x,v,c), 



therefore 



We use (32) and find that 



g„(l — 4:v'^)vg{x, v, c)dv 



V2 



\^(v)rsm{na{v))— + vO{avri 
a{v) 



g{x, V, c)dv — > 0, 



(33) 



fn{x) — c = 2 I g„(l — 4t>^)(yf(x, f , c)df + 2 / g„(l — 4t>^)t>5f(x, f , c)df + o(l), n — +oo. (34) 



as 72 — +00, due to Lemma 2(iii) (note that the function v H- g{x,v,c) is integrable and the function 
v/a{v) is bounded). This fact combined with formula (34) give us (29). □ 



Proof of Theorem 1, (i) and (ii): Theorem 1, (i) follows directly from Proposition 4. Let us prove 
Theorem 1, (ii). Formulas (29) and (33) give us 

e 

^ / X ^ /"i^/ M „ sin(?T,a(f )) , , , 

fn{x)-c = 2 / - ^ ^^ g{x,v,c)dv + o{l) 



e 

= 2 1 \av)r 
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sin(na(t>)) x 



V g{x,v,c) 



a(vi 



dv + o(l). 



The functions sin{na{v)) and v/a{v) are bounded and the function v h- )■ g(x,v,c)/v is integrable over 
[0, e), thus applying Lemma 2(iii) we see that the integral in the right-hand side of the above formula 
converges to zero as n ^ +00, therefore fn{x) — )■ c as n — )■ +00. □ 

Proof of Theorem l,(iii): Assume that for some 6 > the function / has bounded variation on the 
interval [x — 6,x + 6]. We will need the following simple property of functions of bounded variation: 
U f{y) has bounded variation for y G [a,b] and g{x) is monotone for x G [c,d\ and satisfies g{[c,d\) C 
[a,b], then f{g{x)) has bounded variation for x G [c,d]. 

Using the above property we see that there exists ei G (0, 1/4) such that both functions v f{—x log2(^ — 
v)) and V /(— xlog2(| + v)) have bounded variation in the interval v G [0, ei]. 

Formula (28) shows that there exists €2 G (0,1/4), such that a"{v) > for all v G (0,62). We set 
€3 = min(ei, €2) and take e < £3 to be a small positive number (to be specified later). We rewrite formula 
(29) where e is chosen in the above way and c = ^ (/(x + 0) + f{x — 0)): 

e 

fn{x)-c = J le(^)r" ''''|^'|"|''^^ [/(-xlog2(^-^))-/(x + 0)] d.; (35) 


+ / le(^)|-" ''''^7!''^^ [f{-x\og,C2+v))-f{x-0)] dv + o{l), n^+oo. 



Since the two functions v 1— )■ /(— a;log2(^ — v)) and v /(— a;log2(| + v)) have bounded variation 
for V G [0,63], there exist four increasing functions hi{v), 1 < i < 4, satisfying hi{0+) = 0, such that 
/(-a;log2(| - v)) - fix + 0) = h,{v) - h2{v) and f{-x \og^{\ + ^)) - /(^ - 0) = h^{v) - h^{v) for all 
V G [0,63]. We rewrite (35) in the form 

fn{x) - C= Jn^l~ Jn,2 + Jn,3- JnA + W +00, (36) 

where 

/"i^/ M r,sin(na(v)) , , 
J a{v) 



Since hi is positive and increasing, we apply the mean value theorem and conclude that there exists 
6i G (0, e) such that 

T f \ft M-nSin(na(t;)) u t \ f \ft M-nSin(na(t;)) , , 

y «(t;) J a{v) 

01 

Recall that a"(0) > for v G [0,e3] and a'(0) = 2^2 (see (28)), which shows that a'{v) is positive and 
increasing for v G [0, e^], which in turn implies that the function |^(f )|~"'/Q;'(f ) is positive and decreasing 
for V G [0,63]. Applying the mean- value theorem again, we find that there exists some 62 G (^i,e) such 



14 



that 



a[v\ 









f 


X 









01 



sm(na(v)) , 
r-. a (v) 



a(vi 



av = — , / — a [v)dv 



01 



na{92) 



u 



na{6i) 



(39) 



where Si{y) := sm{u)u ^du is the sine-integral function. Since 1^(6*1)1 > 1 and a'{9i) > 2\/2 and 
|Si(x)| < vr for all x > 0, we conclude from (37), (38) and (39) that 



\Jn,i\ < 

In the same way we can obtain corresponding estimates for J„^j, i = 2,3,4, therefore (36) gives us 



vr 



limsup \fn{x) - c| < -^{hi{e) + /12(e) + /13(e) + ^{e)). 



n—^+oo 



Sine hi{0+) = 0, the right-hand side can be made arbitrarily small provided that e is small enough. This 

shows that /„(x) — )■ c as n — )■ +00. □ 
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